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An  Inexpensive  Model  System  for  in  situ  Evaluations  of  Ecoiogicai  Risk  in 
Petrochemical-Contaminated  Terrestrial  Environments 

EXECUTIVE  SUMMARY 

The  overall  objective  of  this  research  was  further  elucidate  the  mechanisms  by 
which  petrochemical  contaminants  influence  terrestrial  ecosystems  and  to  develop 
reliable  in  situ  biomonitoring  tools  for  assessing  risks  to  mammalian  systems  and 
ecosystems.  Specifically,  we  approached  this  objective  by  examining  the  feasibility  of 
using  internal  parasite  communities  harbored  by  resident  small  mammals  as  an  in  situ 
assay  system  for  evaluating  impacts  of  petrochemical  contamination  on  complex 
ecological  processes  in  terrestrial  environments.  Host  parasite  community  ecology 
represents  a  culmination  of  a  myriad  complex  ecological  processes,  including 
interactions  among  hosts,  intermediate  hosts,  and  parasites,  host  immunity,  direct 
toxicity  to  parasites,  and  other  processes. 

Overall  effects  of  petrochemical  soil  contamination  on  the  structure  (species 
composition  and  abundances)  of  the  parasitic  communities  of  cotton  rats  were  evident 
in  the  single-linkage  cluster  analysis  of  similarity  indices.  Host  animals  inhabiting  the 
reference  sites  harbored  a  richer  parasitic  fauna  than  hosts  residing  on  the 
contaminated  sites,  suggesting  that  contamination  negatively  impacted  community 
structure.  Higher  relative  densities  of  surface-dwelling  macroarthropod  groups  were 
observed  on  the  contaminated  sites.  The  most  notable  difference  was  the  terrestrial 
isopods,  which  had  densities  180  times  higher  on  contaminated  than  reference  sites. 
Significant  differences  were  also  noted  in  macroarthropod  diversity,  and  the  results  of 
similarity  comparisons  of  macroarthropod  communities  residing  on  the  contaminated 
and  reference  sites  indicated  community  abundances  and  composition  on  the 
contaminated  sites  were  more  similar  to  one  another  than  to  the  communities  residing 
on  the  reference  sites.  The  results  of  our  survey  suggest  that  helminth  parasite 
communities  of  cotton  rats  and  macroarthropod  communities  in  the  habitat  could  be 
valuable  ecosystem-level  bioindicators  for  evaluating  the  effects  of  soil  contamination 
on  terrestrial  systems. 
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RESEARCH  OBJECTIVES 


Objective  I 

To  evaluate  the  effects  and  sensitivity  of  various  population  and  community 
measures  of  organization  for  parasitic  communities  harbored  by  resident  small- 
mammalian  hosts  inhabiting  terrestrial  environments  contaminated  with  compiex 
mixtures  of  petrochemicals. 

a)  Monitor  seasonal  (winter,  summer)  changes  in  the  prevalence,  over¬ 
dispersion,  abundance,  and  intensity  of  infection  for  each  species  of  parasite  infecting 
resident  small  mammals  collected  from  contaminated  and  uncontaminated  reference 
sites. 

b)  Describe  seasonal  changes  in  organization  (diversity,  species  richness)  and 
compare  structural  similarities  of  parasite  communities  among  ail  study  sites. 

Objective  11 

To  examine  how  abundance  and  diversity  of  host  communities  and  availabiiity  of 
intermediate  parasite  hosts  influence  the  seasonal  diversity,  structure,  and  composition 
of  parasite  communities  on  contaminated  environments. 

There  has  been  no  change  in  the  research  objectives  listed  in  the  original  proposal. 


RESEARCH  APPROACH 


Study  Area 

This  study  was  conducted  at  the  former  Oklahoma  Refining  Company  site,  which 
is  located  in  the  city  of  Cyril,  Oklahoma.  The  site  has  been  an  Environmental  Protection 
Agency  Superfund  waste  site  since  the  refinery  filed  bankruptcy  and  closed  in  1984. 
Three  contaminated  sites  were  utilized;  storage  pit  site,  land  adjacent  to  several 
unlined  asphalt  storage  pits;  oil  sludge  trap  site,  land  adjacent  to  a  series  of  oil  sludge 
sedimentation  ponds;  and  soil  farm  site,  a  3.4  ha  land  farming  site,  which  was  used  to 
process  oil  sludge  waste  materials.  Two  reference  areas  were  also  utilized:  site  one, 
which  was  located  on  property  which  was  owned  by  the  Oklahoma  Refining  Company 
but  not  utilized  in  refining  operations;  and  site  two,  located  on  private  property.  Both 
areas  were  chosen  based  on  their  proximity  to  the  contaminated  sites  (both  areas  were 
located  within  7.2  km  of  the  refinery)  and  their  ecological  similarity.  All  sites  were 
dominated  by  disturbance-adapted  species  such  as  johnsongrass  {Sorghum 
halapense),  Sumac  {Rhus  spp.),  several  species  of  brome  grasses  {Bromus  spp.), 
ragweed  {Ambrosia  spp.),  and  sagewort  {Artemisia  spp.).  Although  cotton  rats  were  by 
far  the  most  commonly  trapped  species  of  small  mammal  on  all  sites,  we  commonly 
trapped  several  other  small  mammal  species,  including  Peromyscus  maniculatus, 
Peromyscus  leucopus,  Reithrodontomys  fulvescens,  and  Mus  musculus. 

Surveys  of  the  major  soil  contaminants  on  the  site  have  found  heavy  metals  and 
organic  waste  materials  from  the  refining  process  on  all  contaminated  sites.  Organics 
were  found  only  in  soils  collected  from  contaminated  sites,  and  the  primary 
contaminants  were  toluene,  xylenes,  pyrene,  anthracene,  napthalene,  phenanthrene, 
benzo(g,h,i)  perylene  and  benzo(a)  pyrene. (Stanley  Engineering  Inc.  1985;  USEPA 
Proj.  No.  W68439).  Heavy  metals  were  also  found  in  soil  samples  from  all 
contaminated  areas.  Lead  and  Chromium  were  detected  from  all  three  areas,  with 
levels  ranging  from  24  to  2700  ppm  for  chromium  and  14  to  304  ppm  for  lead.  Arsenic 
was  detected  in  soil  samples  from  the  storage  pit  and  oil  sludge  trap  sites,  with 
averages  of  104  ppm  and  3  ppm,  respectively.  Aluminum,  barium,  zinc,  and  mercury 
were  only  found  in  soil  samples  from  storage  pits  site,  with  levels  ranging  from  19  to 
53,800ppm  (Stanley  Engineering  Inc.  1985;  USEPA  Proj.  No.  W68439). 

Data  Collection 

In  order  to  examine  possible  differences  in  seasonal  fluctuations  in  the  helminth 
communities  of  cotton  rats  between  contaminated  and  reference  sites,  animals  were 
collected  during  five  time  periods  (September  1993,  February  1994,  August  1995,  April 
1996,  and  October  1996).  Cotton  rats  were  collected  by  removal  trapping  with  snap 
traps  baited  with  peanut  butter.  Animals  were  subsequently  weighed,  sexed, 
necropsied,  and  small  intestinal  tracts  were  removed  and  preserved  in  70%  ethanol  or 
10%  formalin.  Lateral  incisions  were  made  along  the  intestinal  tracts  for  gross 
examinations  and  cestodes  that  were  found  were  removed.  We  removed  the  intestinal 
contents  by  scraping  the  mucosa  with  a  glass  slide  and  examined  the  contents  with  a 
dissecting  microscope  to  find  the  smaller  nematodes.  Parasites  were  examined  and 


taxonomically  identified  using  lactophenol  wet  mounts  for  nematodes  and  borax- 
carmine  stain  mounts  for  cestodes;  total  recovery  of  all  helminth  species  was 
attempted. 

We  sampled  both  macroarthropod  and  microarthropod  communities  during  late 
Summer  and  Fall  1995.  Sampling  of  the  macroarthropod  communities  was  replicated 
three  times  and  we  a  3x4  grid  of  pit-fall  traps  on  each  site.  All  arthropods  were 
removed  at  the  end  of  the  sampling  were  identified  and  enumeration  to  family,  and 
categorized  into  general  taxonomic  groups  for  relative  density  comparisons  between 
the  contaminated  and  reference  sites.  Other  macroarthropod  groups  were  identified 
into  Classes  (Isopoda,  Symphyla,  Chilopoda,  and  Diplopoda).  Sampling  of  the 
microarthropod  communities  was  replicated  5  times  from  late  August  to  early  October 
1995.  During  each  sampling  period  we  removed  four  soil  cores  from  each  study  site. 
Soil  cores  were  immediately  wrapped  in  aluminum  foil  and  stored  on  ice  to  reduce 
arthropod  movements  until  they  could  be  returned  to  the  laboratory  within  2-4  hours 
(Seastedt  and  Crossley  1981)  for  Tullgren  extraction  (Merchant  and  Crossley  1970)  for 
7  days.  Microarthropods  were  collected  and  mites  (Acarina)  were  identified  to 
Suborder,  and  other  microarthropods  were  identified  to  Order. 

Data  Analysis 

We  used  a  completely  randomized  design  with  repeated  measures  for 
comparisons  of  abundance,  intensity  and  overdispersion.  The  comparisons  were  made 
using  PROC  MIXED  (SAS,  1996)  with  sources  of  variation  including  treatment,  site 
within  treatment  (error  term  for  treatment),  sampling  date,  treatment  by  date  interaction, 
and  the  residual.  A  compound  symmetric  model  was  used  to  model  the  covariance 
structure  of  the  repeated  measurements.  If  the  treatment  by  sampling  date  interaction 
was  significant,  simple  effects  of  treatment  were  analyzed  using  the  SLICE  option  for 
the  LSMEANS  statement.  Satterthwait’s  approximation  was  used  for  calculation  of  the 
degrees  of  freedom  of  the  pooled  error  term.  If  the  treatment  by  sampling  date  was  not 
significant,  the  main  effects  were  analyzed  with  the  DIFF  option.  Comparisons  of 
prevalence  were  made  utilizing  Fisher’s  Exact  Test  for  detection  of  heterogeneity 
among  the  five  trapping  areas  during  each  season  (PROC  FREQ,  SAS  1996). 

Communities  for  each  host  population  were  described  by  measures  of  diversity, 
mean  species  richness,  and  similarity.  Diversities  were  calculated  by  using  the 
complement  of  Simpson’s  index,  1-  D  (Krebs,  1989).  Comparisons  of  mean  species 
richness  and  species  diversity  among  treatments  were  made  using  analysis  of  variance 
with  repeated  measures  as  described  previously.  Similarities  in  species  composition 
and  relative  abundances  of  helminth  communities  among  cotton  rat  populations  were 
calculated  using  Horn’s  index  (Horn,  1966)  and  relationships  between  communities 
were  depicted  using  single-linkage  cluster  diagrams  (Krebs,  1989).  Statistical 
significance  for  all  hypothesis  tests  was  set  a  priori  at  P  <  0.05. 


RESEARCH  OVERVIEW 


Helminth  Fauna  and  Prevalence 

All  but  four  of  the  340  cotton  rats  which  we  collected  were  infected  with  at  least 
one  species  of  helminth  (Table  1).  Three  species  of  cestodes  {Raillietina  sigmodontis, 
Schizotaenia  sigmodontis,  and  Hymenolepis  dimunuta)  and  four  species  of  nematodes 
(Protospirura  muris,  Longistriata  adunca,  Syphada  sigmodontis,  and  Strongyloides 
sigmodontis)  were  recovered.  Longistriata  adunca,  a  trichostrongylid  nematode,  was 
the  most  prevalent  species,  infecting  316  of  the  340  cotton  rats.  Schizotaenia 
sigmodontis  was  the  most  prevalent  cestode  in  the  study,  infecting  49%  of  the  cotton 
rats  which  we  surveyed. 

Prevalences  among  helminth  populations  within  each  season  were  extremely 
variable,  and  no  consistent  patterns  of  difference  between  the  three  contaminated  and 
two  reference  sites  were  evident  for  most  species  (Figs.  1  and  2).  Longistriata  adunca 
was  found  in  over  80%  of  the  animals  on  all  sites  in  the  first  four  seasons,  but 
prevalence  dropped  below  60%  in  animals  from  the  storage  pit  site  in  October  1995 
(Fig.  2).  Hymenolepis  dimunuta  was  not  observed  in  September  1 993  or  February 
1994,  and  Syphada  sigmodontis  v/as  not  observed  in  October  1995;  all  other  species 
were  observed  in  every  sampling  period.  There  were  three  seasons  in  which  the 
prevalence  of  a  helminth  species  on  both  reference  sites  was  significantly  different  from 
any  of  the  contaminated  sites.  Syphada  sigmodontis  was  more  prevalent  (P  <  0.05)  in 
cotton  rats  collected  from  the  reference  sites  than  the  sludge  trap  and  soil  farm  sites  in 
April  1995.  In  September  1994  prevalence  of  infection  with  Strongyloides  sigmodontis 
was  greater  (P  <  0.05)  in  cotton  rats  from  the  reference  sites  than  the  storage  pit  and 
sludge  trap  sites,  where  this  species  was  not  recovered  from  any  host  animals 
examined.  Protospirura  muris  was  five  times  more  prevalent  on  the  reference  sites  than 
the  sludge  trap  and  soil  farm  sites  (P  <  0.05)  in  October  1995.  Although  consistent 
differences  in  prevalences  from  one  season  to  another  between  contaminated  and 
reference  sites  were  lacking  for  all  species,  the  differences  in  prevalences  of 
Strongyloides  sigmodontis  among  the  contaminated  study  sites  over  the  course  of  the 
survey  were  noteworthy.  Prevalences  of  Strongyloides  sigmodontis  varied  widely  on  the 
reference  sites  across  seasons,  ranging  from  10  to  50%,  while  remaining  consistently 
low  (0  to  15%)  in  cotton  rats  from  all  of  the  contaminated  sites  throughout  the  survey. 

Helminth  Abundance 

Total  abundances  of  cestodes  in  cotton  rats  in  our  survey  were  influenced 
significantly  by  season  of  collection  (P  <  0.0001)  but  not  treatment.  The  differences  in 
cestode  abundances  were  primarily  due  to  low  overall  cestode  numbers  in  February 
1994  and  ^ril  1995  compared  to  the  fall  collections.  The  abundances  of  Schizotaenia 
sigmodontis  infections  were  influenced  by  season  of  collection  and  by  treatment  (Table 
1 ).  Schizotaenia  sigmodontis  was  twice  as  abundant  in  host  animals  on  the  reference 
sites  than  the  contaminated  sites  (P  <  0.005).  The  influence  of  season  of  collection  (P  < 
0.01)  on  the  abundance  of  Schizotaenia  was  primarily  due  to  the  50%  decrease  in 
abundances  of  the  species  during  October  1 995  compared  to  other  seasonal 


collections.  Abundances  of  Raillietina  sigmodontis  infections  were  not  influenced  by 
treatment  (P  >  0.05),  but  were  significantly  influenced  by  season  of  collection  (P  < 
0.0001),  and  there  was  also  a  significant  season  by  treatment  interaction  (P  <  0.05). 
Analysis  of  the  simple  effects  indicated  that  there  was  a  significantly  larger  variation  in 
abundances  Raillietina  sigmodontis  on  the  contaminated  sites,  with  higher  abundances 
in  the  Fall  collections  and  lower  abundances  in  April  1 995  and  February  1 994,  when  no 
infected  animals  were  recovered  from  the  contaminated  sites.  Although  a  similar 
seasonal  pattern  of  higher  fall  abundances  was  present  in  animals  from  the  reference 
populations,  variation  in  abundances  of  Raillietina  sigmodontis  on  the  reference  sites 
was  considerably  lower  and  was  only  significantly  (P  <  0.05)  different  from  other 
seasons  when  it  increased  in  October  1995.  Abundances  of  infection  with  Hymenolepis 
dimunuta  were  significantly  different  across  seasons  (P  <  0.0005).  There  was  also  a 
significant  treatment  by  season  interaction  (P  <  0.05)  due  to  higher  abundances  of 
Hymenolepis  dimunuta  in  cotton  rats  from  the  reference  sites  in  October  1 995. 

Remarkable  differences  were  also  observed  in  seasonal  abundances  of  several 
nematode  species  infecting  cotton  rats  in  our  study.  Abundances  of  Protospirura  muris 
showed  significant  seasonal  fluctuations  (P  <  0.0001),  with  a  strong  treatment  by  date 
interaction  (P  <  0.005).  Seasonal  changes  in  abundances  of  Protospirura  muris 
infections  in  cotton  rats  from  reference  sites  varied  almost  8-fold  across  seasons,  while 
significant  seasonal  changes  in  abundances  on  contaminated  sites  were  not  apparent 
(Table  1).  Similar  to  Protospirura  muris,  abundances  of  Syphacia  sigmodontis  infections 
were  significantly  different  across  seasons  (P  <  0.05) ,  with  a  significant  season  by 
treatment  interaction  (P  <  0.01 ).  Reference  host  populations  had  significantly  greater 
fluctuations  in  abundances  of  Syphacia  sigmodontis  infections  across  seasons 
compared  to  host  populations  from  contaminated  sites.  Abundances  of  Syphacia 
sigmodontis  from  reference  host  populations  were  lower  in  February  1 994  and  April 
1995  than  September  1994  and  October  1995,  while  abundances  of  infections  in  host 
animals  from  contaminated  sites  were  only  different  from  other  seasonal  collections  in 
September  1 994.  The  large  mean  abundance  of  Syphacia  sigmodontis  infections  from 
hosts  on  the  contaminated  sites  in  September  1994  was  due  to  a  few  heavily  infected 
animals,  including  one  host  which  harbored  over  9000  individuals.  Abundances  of 
Strongyloides  sigmodontis  were  also  influenced  by  season  (P  <  0.01 ),  with  a  strong 
treatment  by  season  interaction  (P  <  0.005).  Analysis  of  the  simple  effects  showed  that 
seasonal  abundances  of  Strongyloides  sigmodontis  infections  also  varied  widely  in  host 
populations  on  reference  sites,  peaking  in  September  1 993  and  1 994;  abundances 
remained  very  low  (less  than  0.1  worms  per  host)  in  cotton  rat  populations  from 
contaminated  sites  across  all  seasons. 

Helminth  Community 

Mean  species  richness  of  gastrointestinal  parasite  communities  within  host 
populations  was  significantly  influenced  by  treatment  (P  <  0.05;  Table  2).  Host  animals 
on  the  reference  sites  consistently  supported  a  greater  average  number  of  helminth 
species  than  those  from  contaminated  sites  in  all  seasonal  collections.  Mean  species 
richness  also  differed  across  seasons  (P  <  0.05),  being  significantly  higher  in 
September  1994  than  in  April  and  October  1995.  Diversity  (complement  of  Simpson’s 
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index)  did  not  appear  to  change  seasonally  or  differ  significantly  between  treatments  (P 
>  0,05),  ranging  from  0.232  to  0.603. 

Comparisons  of  the  similarity  of  species  composition  and  abundances  among 
helminth  communities  showed  that  the  reference  sites  were  more  apt  to  cluster  together 
than  with  the  contaminated  sites,  especially  in  February  1994,  September  1994,  and 
April  1995  (Fig.  3).  Although  reference  site  1  was  most  similar  to  the  sludge  pits  site  in 
October  1995,  the  two  reference  sites  clustered  together  at  >  97%  similarity.  In 
September  1993  the  three  contaminated  sites  were  clustered  together,  but  the 
reference  sites  were  not. 

An  three  cestodes  and  one  of  the  nematode  species  recovered  from  cotton  rats 
in  this  study  have  complex  life  cycles  which  require  arthropod  intermediate  hosts  and 
all  four  of  these  species  demonstrated  a  significant  treatment  effect  or  treatment  by 
season  interaction.  It  seems  plausible  that  some  of  this  difference  is  the  result  of 
impacts  (positive  or  negative)  from  soil  contamination  on  the  availability  of  intermediate 
hosts  (see  below  in  Arthropod  section),  which  could  play  a  profound  role  in  the 
regulation  of  these  helminth  species  in  their  hosts. 

The  other  three  nematode  species  we  recovered  all  have  direct  life  cycles, 
although  there  are  still  some  notable  life  cycle  differences  among  them.  Longistriata 
adunca  and  Syphacia  sigmodontis  are  both  shed  by  the  host  directly  into  the 
environment,  where  they  are  transmitted  to  new  hosts  by  ingestion  of  eggs  or  larvae. 
The  results  of  our  study  suggest  that  seasonal  fluctuations  in  abundances  of  infection 
with  Syphacia  sigmodontis  were  potentially  influenced  by  contamination,  but  that 
abundances  of  Longistriata  adunca  were  more  significantly  influenced  by  seasonal 
factors  and  not  by  contamination.  Species  of  the  genus  Strongyloides  have  a  unique 
heterogenic  life  cycle,  with  alternation  of  free-living  and  parasitic  female  generations. 
The  males  are  not  parasitic,  residing  in  a  free-living  form  in  the  soil.  The  parasitic 
filariform  larvae  enter  the  host  by  direct  penetration  of  the  skin  (Melvin  and  Chandler, 

1 950).  It  seems  probable  that  this  life  cycle  would  expose  the  free-living  generations 
directly  to  the  contaminants  in  the  soils,  which  could  have  a  negative  effect  on  the 
densities  of  free-living  generations  and  on  the  availability  of  males,  which  would  explain 
the  consistently  low  prevalences  and  abundances  of  parasitic  females  we  observed  on 
contaminated  areas. 

Overall  effects  of  petrochemical  soil  contamination  on  the  structure  (species 
composition  and  abundances)  of  the  parasitic  communities  were  evident  in  the  single¬ 
linkage  cluster  analysis  of  similarity  indices.  Also,  host  animals  inhabiting  the 
reference  sites  harbored  a  richer  parasitic  fauna  than  hosts  residing  on  the 
contaminated  sites,  suggesting  that  contamination  negatively  impacted  community 
structure.  These  results  suggest  that  soil  contamination  was  an  important  contributor  in 
shaping  helminth  communities  of  the  cotton  rat  host  populations.  Much  of  the 
variability  in  the  degree  of  impact  of  contamination  was  probably  due  to  differences  in 
the  life  cycles  of  the  individual  parasite  species.  Given  the  differences  in  populations 
and  overall  community  structure  which  we  observed,  host-parasite  relationships  seem 
to  have  promise  as  a  biomonitoring  tool  for  the  assessment  of  community-level  effects 
of  contamination  in  terrestrial  ecosystems. 


Arthropod  Community 

Relative  densities  of  the  total  macroarthropod  communities  were  greater  (P  < 
0.005)  on  the  contaminated  than  reference  sites  and  aiso  varied  significantly  (P  < 
0.0001)  across  the  three  collection  dates  (Tabie  3).  A  significant  treatment  by  sampling 
date  interaction  (P  <  0.0001 )  was  indicated  because  total  macroarthropod  relative 
densities  decreased  62  to  71%  on  contaminated  study  sites  by  the  24  October  survey 
compared  to  the  two  earlier  surveys,  while  the  relative  densities  on  the  reference  sites 
remained  stable.  Relative  densities  of  hymenopterans  were  50  to  100  %  greater  (P  < 
0.01 )  on  contaminated  than  reference  sites  in  the  October  surveys.  Relative  densities 
of  coleopterans  (beetles)  were  not  different  (P  >  0.05)  between  sampling  dates  or 
treatments  remaining  fairly  constant  across  all  surveys  (Table  3).  Although  the  main 
effects  of  sampling  date  and  treatment  did  not  (P  >  0.05)  affect  relative  densities  of 
total  arachnids,  there  was  a  significant  treatment  by  sampling  date  interaction  (P  < 
0.0005).  Relative  densities  of  arachnids  were  four  fold  greater  on  contaminated  sites 
compared  to  reference  sites  on  24  October.  Relative  densities  of  Order  Araneae 
showed  a  significant  (P  <  0.05)  treatment  by  sampling  date  interaction  due  to  lower 
densities  on  the  contaminated  sites  on  1 0  October  and  on  the  reference  sites  on  24 
October.  In  both  of  these  surveys  the  densities  of  Order  Araneae  were  50  to  60% 
lower  on  these  sites  than  the  other  surveys  of  the  reference  and  contaminated  sites. 

Overall  isopod  relative  densities  were  1 80-fold  greater  (P  <0.005,  Fig.  4)  on  the 
contaminated  sites  than  on  the  reference  sites.  A  significant  treatment  by  sampling  date 
interaction  (P  <  0.0001 )  was  also  indicated  as  the  relative  densities  of  isopods  on  the 
contaminated  sites  decreased  over  80%  on  24  October  compared  to  the  29  August  and 
10  October  surveys.  Relative  densities  of  isopods  remained  consistently  low  (an 
average  of  less  than  one  isopod  per  trap)  on  the  reference  sites  throughout  all  survey 
periods. 

Comparisons  of  the  densities  of  eight  microarthropod  groups  showed  that  soil 
contaminants  had  little  or  no  measurable  impact  on  the  densities  of  populations  in  this 
study  (Table  4) 

Comparisons  similarity  and  diversity  of  macroarthropod  communities  indicated 
that  soil  pollution  had  significant  effects  on  community  structure.  Comparisons  of 
relative  abundances  and  species  composition  of  macroarthropod  communities  on 
contaminated  and  reference  sites  revealed  that  the  contaminated  sites  clustered 
together  for  all  three  of  the  collections,  with  the  discrimination  between  the  reference 
and  contaminated  sites  most  pronounced  on  29  August  and  10  October  (Fig.  5).  In  the 
29  August  and  10  October  surveys  the  contaminated  sites  all  clustered  together  at  > 
90%  similarity,  while  similarity  between  the  reference  and  contaminated  sites  was 
beiow  50%. 

Our  results  suggest  that  macroarthropod  communities  are  considerably  more 
sensitive  to  the  complex  mixtures  of  contaminants  associated  with  oil  refinery  waste 
than  microarthropod  communities.  The  most  remarkable  difference  that  we  observed 
was  undoubtedly  the  enormous  differences  in  the  densities  of  terrestrial  isopods 
between  contaminated  and  reference  areas.  Isopods  were  the  dominant  group  on  the 
contaminated  sites,  frequently  comprising  >  90%  of  the  total  arthropods  surveyed. 
Although  isopods  occurred  on  the  reference  sites,  there  were  never  more  than  two 
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individuals  in  a  sample.  This  enormous  difference  in  relative  densities  could  have  been 
caused  by  exposure  to  certain  toxicants  on  the  site,  some  of  which  have  been 
demonstrated  to  increase  reproductive  allocation  and  output  in  certain  isopod  species. 
For  example,  Van  Brummelen  et  al.  (1996a)  observed  a  significant  stimulation  of 
reproduction  in  laboratory  populations  of  Oniscus  ase//us  following  exposure  to 
fluoranthene,  phenanthrene,  benz[a]anthracene,  and  ben20[a]pyrene,  which  isopods 
do  accumulate  in  the  environment  (Van  Brummelen  et  al.  1996b).  Donker  et  al.  (1993) 
observed  early  onset  of  sexual  maturation  and  increased  reproductive  allocation  in 
populations  of  the  isopod  Porcellio  scaber  which  inhabited  areas  near  a  zinc  smelter 
and  a  lead  mine  compared  to  populations  in  uncontaminated  areas.  Porcellio  scaber 
populations  have  also  been  shown  to  adapt  to  high  levels  of  cadmium,  which  was  also 
observed  to  stimulate  growth  in  isopod  populations  which  resided  in  close  proximity  to 
a  zinc  smelter  (Donker  and  Bogert  1991).  In  contrast.  Hunter  et  al.  (1987)  observed 
lower  population  densities  of  isopod  populations  in  grassland  areas  contaminated  with 
copper  and  cadmium.  It  seems  likely  that  the  high  densities  of  isopods  we  observed  on 
the  contaminated  sites  in  our  study  were  due  to  increased  recruitment  rates,  although 
indirect  effects  of  contamination  (such  as  direct  effects  of  contamination  on  predators  of 
isopods)  could  not  be  ruled  out  as  well. 
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Table  2,  Indices  of  species  richness  (mean  number  of  helminth  species  per  host  animal)  and  community 
diversity  of  gastrointestinal  helminth  communities  of  hispid  cotton  rats  inhabiting  the  Oklahoma  Refining 
Co.  Super fund  waste  site^  Caddo  County^  Oklahoma,  from  September  1993  to  October  1995. 
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Table  3.  Mean  relative  population  densities  (average  number  of  individuals  per  trap)  of  major  macroarthropod  group; 
collected  from  three  contaminated  and  two  reference  sites  on  the  Oklahoma  Refining  Co.  Superfund  waste  site,  Caddo 
County,  Oklahoma,  August  to  October  1995. 
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ble  4-  Mean  population  densities  (individuals  x  10^/m^)  of  major  microarthropod  groups  collected  from  three 
ntaminated  and  two  reference  areas  of  the  Oklahoma  Refining  Co.  Superfund  waste  site,  Caddo  County,  Oklahoma, 
lues  represent  pooled  data  across  surveys  from  August  to  October  1995. 
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Legend 


Figure  2:  Seasonal  prevalences  of  nematode  species  collected  from  cotton  rats  on  the 
three  contaminated  and  two  reference  sites  on  the  Oklahoma  Refining  Co.  Superfund 
Waste  Site,  Caddo  Co.,  Oklahoma,  August  1993  to  October  1995.  Different  letters 
above  the  bars  denote  significant  differences  in  prevalence  among  study  populations 
within  each  season  (calculated  with  Fisher’s  Exact  Test,  P  ^  0.05). 
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Figure  3:  Single-linkage  duster  diagrams  depicting  the  similarity  (Horn  1966)  of 
helminth  communities  surveyed  on  three  contaminated  and  two  reference  sites  on  the 
Oklahoma  Refining  Company  superfund  waste  site,  Caddo  County,  Oklahoma,  August 
to  October  1995. 
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Figure  4:  Relative  densities  (average  number  of  individuals  per  trap)  of  terrestrial 
isopods  recovered  from  contaminated  and  reference  sites  on  the  Oklahoma  Refining 
Company  superfund  site,  Caddo  County,  Oklahoma,  August  to  October  1995. 


Figure  5:  Single-linkage  cluster  diagrams  depicting  the  similarity  (Horn  1966)  of 
macroarthropod  communities  surveyed  on  three  contaminated  and  two  reference  sites 
on  the  Oklahoma  Refining  Company  superfund  waste  site,  Caddo  County,  Oklahoma, 
August  to  October  1995. 
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